Thickness dependence of linear and quadratic magneto-optical Kerr effects in 

ultrathin Fe(OOl) films 



OO 
O 

o 

(N 



c/5 



a 

c 

o 

(N 
> 

ON 
On 

in 

cn 
o 

oo 
O 



x 



M. Buchmeier^ R. Schreiber, D. E. Biirgler, and C. M. Schneider 
Institut fiir Festkorperforschung (IFF-9) and JARA-FIT, 
Forschungszentrum Jiilich GmbH, D-52425 Jiilich, Germany 
(Dated: May 7, 2008) 

Magneto-optical Kerr effect (MOKE) magnetometry is one of the most widely employed techniques 
for the characterization of ferromagnetic thin film samples. Some information, such as the magnitude 
of coercive fields or anisotropy strengths can be readily obtained without any knowledge of the 
optical and magneto-optical (MO) properties of the material. On the other hand, a quantitative 
analysis, which requires a precise knowledge of the material's index of refraction n and the MO 
coupling constants K and G is often desirable, for instance, for the comparison of samples, which 
are different with respect to ferromagnetic layer thicknesses, substrates, or capping layers. While 
the values of the parameters n and the linear MO coupling parameter K reported by different 
authors usually vary considerably, the relevant quadratic MO coupling parameters G even for Fe are 
completely unknown. Here, we report on measurements of the thickness dependence (0-60 nm) of the 
linear and quadratic magnetooptical effects in epitaxial bcc-Fe(OOl) wedge-type samples performed 
at a commonly used laser wavelength of 670 nm. By fitting the thickness dependence we are able to 
extract a complete set of parameters n, K, (Gn — G12), and G44 for the quantitative description 
of the MOKE response of bcc-Fe(OOl). We find the parameters n, K, and G to significantly differ 
for films thinner than about 10 nm as compared to thicker films, which is indicative of a thickness 
dependence of the electronic properties or of surface contributions to the MOKE. The magnitude 
of the quadratic magnetooptical effect is found to be about one third of the record values reported 
recently for Co2FeSi. 

PACS numbers: 78.20.Ls 78.20.Ci 78.20.Bh 75.70.Ak 75.30.Gw 75.60.-d 



I. INTRODUCTION 

The magneto-optical Kerr effect (MOKE) of ferromag- 
netic (FM) thin films has been a field of intensive stud- 
ies over the last three decades. This interest in MOKE 
was motivated by three aspects: (i) its importance as 
a experimental magnetometric tool, (ii) as a means to 
measure the band structure of FM materials, and (hi) 
its application in magneto-optical (MO) storage media. 
MOKE is probably the tool most widely employed for 
the magnetometric characterization of thin-film samples 
employed for spintronics. Among its most common ap- 
plications are the quantitative determination of the coer- 
civity, magnetic anisotropy, and interlayer exchange cou- 
pling from the analysis of hysteresis loops recorded with 
the MOKE signal. Other prominent applications are the 
investigation of spin dynamics in the time-domain and 
magnetic domain imaging. The main advantages of the 
MOKE over other techniques are its compatibility with 
high magnetic fields, surface sensitivity with a typical in- 
formation depth of some 10 nm, a time resolution down 
to the sub-picosecond regime, a reasonable spatial reso- 
lution of the order of about 0.5 /im, and robust and in- 
expensive experimental setups. 

Many applications neglect the absolute magnitude of 
the Kerr effect, which is given by the magnitude and 
phase of the complex Kerr angle. Instead they describe 
the dependence of the normalized Kerr angle on the 
magnetization angle by means of adjustable response co- 
efficients. This type of description has the advantage 
that it does not require any knowledge of the materi- 



al's optics, yet it is sufficient to extract a lot of infor- 
mation, such as the magnitude of the magnetocrystallinc 
anisotropy of FM single-layer systems, 1 :- or the antifer- 
romagnetic interlayer exchange coupling of FM double- 
layer systemsi^ On the other hand, the absolute magni- 
tude of the complex Kerr angle provides valuable infor- 
mation, which can support the experimental data, and 
can be employed to quantitatively compare samples with, 
for instance, different FM layer thicknesses, substrates, 
or capping layers and to determine the thickness of the 
FM layers, the MOKE information depth, etc. However, 
the full quantitative MOKE information is generally not 
linked by simple analytic formulae to the material prop- 
erties, which are the indices of refraction n and the lin- 
ear and quadratic MO coupling parameters K and G of 
all involved layers. Even in ultrathin films of only some 
nanometers thickness, a linear dependence of the size of 
the MOKE response on the FM layer thicknesses, known 
as additivity law - which has been claimed by Qiu et al. 
in Ref.d - is generally not validi^ 

Therefore, a general numerical treatment of the MOKE 
by solving Maxwell's equations and the standard bound- 
ary conditions is indispensable for the quantitative in- 
terpretation of the Kerr angle. A prerequisite for this 
calculation is the precise knowledge of the optical and 
magnetooptical material parameters. Although spcctro- 
scopically determined values for n and K are available for 
many materials, the overall agreement of the data from 
different sources is often, as for instance in the case of 
bcc-Fe^ rather poor. The strong variation of the tabu- 
lated optical constants is frequently attributed to surface 
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contamination or oxidation of the ex-situ measured sam- 
ples, but thickness or quality variations of the films are 
also plausible. For instance, (i) in Ref.dthe Kerr angle of 
Fe has been found to oscillate as a function of the thick- 
ness of a capping Au layer due to quantum well states, 
(ii) spectroscopic MOKE data of thin Fe films show fea- 
tures, which cannot be explained by the bulk electronic 
band structure of Fe^ and (iii) there is strong evidence 
for interfacial contributions in the MOKEiiS 

Moreover, the literature values are almost exclusively 
limited to the first order linear MO coupling parame- 
ter K. For many FM materials, e.g. from elementary 
Fe^ i 11 ' 12 ! 13 to the more complex Heusler alloy Co2FeSi^, 
however, the second order quadratic coupling constants 
G have a comparable impact on the Kerr angle. To our 
knowledge the only published value of a second order 
MO coupling constant is the imaginary part of G44/K 
of Ni of about Im{G^j K) = —0.02 at the wavelength 
A = 514 nm, which has been experimentally determined 
from Brillouin light scattering data by Giovannini et al*£> 
The second order magnetooptical parameters, which are 
G\\, G12, and G44 for systems with cubic symmetry, give 
rise to MO effects quadratic and even in the magnetiza- 
tion M. These effects are known as quadratic MOKE 
(QMOKE) or Voigt effect in reflection, and have recently 
received a lot of attention^ ' 11 ' 12 ] 13 ! 14 Effects quadratic 
in M also turn out to be important in magnetization- 
dependent second harmonic generation (MSHG)f 1 ^ x- 
ray magnetic linear dichroism and the closely related 
x-ray Voigt effect^ All the more it is surprising, that 
to our knowledge no one has yet determined a full set 
of corresponding material parameters G to describe the 
QMOKE. 

Apart from their practical significance for MOKE mag- 
netometry, the MO coupling parameters are also impor- 
tant in fundamental research. From a microscopic point 
of view the MO coupling is due to the interplay of the 
exchange exchange interaction leading to a splitting of 
the bands and the spin-orbit (SO) coupling. It is there- 
fore closely related to the magnetocrystalline anisotropy 
energy, which also arises due to SO coupling. The MO 
coupling parameters may be seen as an important probe 
for the fundamental electronic interactions in FM mate- 
rials, e.g. spectroscopic MOKE is a widely used stan- 
dard tool to evaluate the band structure in FM metals. 7 
While the linear MO effect is due to first order SO cou- 
pling, the quadratic MO coupling is thought to be caused 
by second order SO coupling terms. The SO coupling 
is known to be altered at the interfaces and in ultra- 
thin films, which gives rise to the well-known thickness 
dependence of the magnetocrystalline anisotropy energy, 
e.g. in Fe^ and interfacial MOKE contributions, which 
have been observed by various authors* 1 ^ 

In this contribution we report on a magnetometric 
study of the magnetooptical response of bcc-Fe(OOl) 
wedge-type samples with thicknesses ranging from to 
60 nm. We have determined both components of the com- 
plex Kerr angle, the Kerr rotation and the Kerr ellipticity. 



Effects linear and quadratic in M, LMOKE and QMOKE 
respectively, are separated by fitting the hysteresis loops 
to a single domain model. The QMOKE, which is known 
to be anisotropic, i.e. it depends on the sample orien- 
tation with respect to the plane of incidence, has been 
determined for both Fe(001)[110] and Fe(001)[100] direc- 
tions parallel to the plane of incidence. By fitting the 
thickness dependence of LMOKE and QMOKE we are, 
for the first time, able to extract a full set of Fe mate- 
rial parameters n, K, (Gn — G12), and G44 at a light 
wavelength of 670 nm. We find a sizeable thickness de- 
pendence, which, however, seems to be mainly of non- 
magnetic origin. The main effect is an increased index 
of refraction for Fe film thicknesses below about 10 nm 
as compared to thicker films. A large maximum absolute 
value of the quadratic Kerr effect (QMOKE) of 0.37 mrad 
is found at about 22 nm Fe thickness. This value is of a 
comparable order of magnitude as the recently reported 
record QMOKE values for Co 2 FeSi^i 

The paper is organized as follows: In Sec. Ill Al we 
describe briefly the sample preparation. Details of the 
experimental MOKE setup and the data recording are 
given in Sec. Ill Bl Scction HiTl deals with the modeling of 
the MOKE and the hysteresis loops. The results are pre- 
sented and discussed in Sec.lIVI Finally, we summarize 
our results in Sec.lVl 



II. EXPERIMENTAL DETAILS 

A. Sample preparation and experimental MOKE 
setup 

Epitaxial Fe(wedge)/Ag(lnm)/Au(2nm) have been 
prepared by molecular beam epitaxy on top of a 
GaAs/Ag(001) buffer system. The Au capping layer has 
has been chosen thick enough to prevent oxidation and 
thin enough to be able to determine large Kerr angles. 
The Ag interface layer has been introduced in order to 
prevent a possible alloying of Fe and Au. The prepara- 
tion is described in detail elsewhere. 19 ' 20 All thicknesses 
have been precisely determined using a calibrated quartz 
crystal monitor. The Fe thickness has been varied con- 
tinuously between and 8 nm for sample A and stepwise 
in sample B with discrete Fe thicknesses of 5, 8, 12, 18, 
24, 32, 44, and 60 nm. 

B. MOKE setup 

The MOKE measurements were performed using light 
from a Toshiba TOLD9231M multi-mode laser diode 
with a wavelength of 670 nm and a spectral half- width of 
less than 2 nm. An in-plane magnetic field with a maxi- 
mum strength of 0.7 T is generated by a Broker electro- 
magnet with FM yoke and measured with a Hall probe. 
The sample is mounted on a manually rotatable holder 
with an angle scale, which allows a determination of the 
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angle with respect to the plane of incidence with a pre- 
cision of about ±2°. 

The optical setup consists of a light beam with its plane 
of incidence parallel to the external field direction illu- 
minating the sample under an angle of incidence of 15° 
with respect to the sample normal (see Fig. [I). The in- 
cident light is polarized in p direction (electric field com- 
ponent in the plane of incidence) employing a rotatable 
Glan- Taylor type polarizer. The phase of the reflected 
light is modulated at a frequency / = 50 kHz with a 
HINDS PEM-90 photo elastic modulator (PEM) at diag- 
onal modulation axis (rotated by 45° with respect to the 
p-direction) and at a retardation of 90° . A quarter wave 
plate with its retardation axis parallel to the p direction 
can be placed between sample and PEM. Consecutively 
the light passes through an analyzer of the same type 
as the polarizer, but oriented in s direction. The light 
intensity is then converted into an electrical voltage by 
a homemade diode detector. In order to avoid possible 
multiple reflections between the sample and the optical 
elements in the reflected beam-path, which might impair 
the data by parasite signals, a diaphragm with a diameter 
of about 2 mm is placed right after the sample. Analyzer, 
PEM, and quarter wave plate are slightly tilted with re- 
spect to the optical axis so that the light reflected back 
to the sample is blocked by the diaphragm. 

With this setup the small / (50 kHz) component deter- 
mined with a Lock-In amplifier is to first order propor- 
tional to the p-ellipticity e times the reflected intensity, 
while the much larger 2/ (100 kHz) and DC components, 
recorded simultaneously with a multimeter are to first 
order proportional to the reflected intensity alone^ By 
introducing a quarter wave plate between sample and 
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FIG. 1: Experimental MOKE setup. The quarter wave plate 
is introduced to measure the Kerr rotation 9 instead of the 
Kerr ellipticity e. For clarity, the diaphragm is not shown. 



PEM we are able to measure the p-rotation 9 instead of 
the ellipticity. The two Kerr angle components (9 and e) 
are calculated by dividing the measured / component by 
the 2/ component. 

As the amplification factor of the detector is frequency- 
dependent, the measured Kerr angle has to be calibrated. 
For this purpose, we have used a thick Au sample and 
turned the analyzer out of the p direction by an angle 
while recording the ellipticity and rotation signal. The 
data are then adjusted such that the measured depen- 
dence of the resulting rotation 9 and ellipticity e as a 
function of the polarizer angle <3> matches the theoretical 
relationship 



0(tf) + ie(tf) = — sin(#), 

Tpp 



(1) 



where r ss and r pp are the diagonal Fresnel reflection co- 
efficients. A value r ss /r pp = 1.0003 — 0.0379z has been 
calculated using the Fresnel formulas and taking the lit- 
erature valued n = 0.100 + 3.653i for the index of refrac- 
tion of Au at our laser wavelength. We have performed 
the calibration procedure for both the real and imagi- 
nary parts of Eq. ifT]). i.e. with and without quarter wave 
plate, and find an excellent agreement for the calibration 
factors (less than 2% difference). This corroborates that 
the calibration works properly and possible detriments of 
the measurement by multiple reflections are indeed well 
under control. 



III. NUMERICAL MODELLING OF THE MOKE 
AND HYSTERESIS LOOPS 



The optical and MO material properties can be de- 
scribed by the permittivity tensor e+j, which can be 
expanded in a power series of the Cartesian direction 
cosines of the magnetization M: 



(M=0) 



G 



tjks m k m s 



(2) 



The number of independent linear and quadratic MO 
coupling constants, Kijk and Gijks respectively, is re- 
duced by the symmetry of the crystal and the On- 
sager principle eij(M) — €ij{—M), For cubic symme- 
try the permittivity tensor is completely defined by five 
quantities: The nonmagnetic part of the permittivity 
e (iW=o) _ n 2^ which j s given by the square of the index of 
refraction n, the linear MO coupling constant K = Kijk, 
and three independent quadratic MO coupling param- 
eters Gnu = Gn, Gujj = G12 and Gyij = G44. In- 
stead of the linear MO coupling K, the Voigt parameter 
Q = iKm/ 'e( M=0 ) is frequently used. The complex Kerr 
angle, which is a measure of the magnitude of the MOKE 
signal can be calculated using the standard 4x4 ma- 
trix formalism as explained in Refs j23ll 2411251 . Our open- 
source computer program developed for the calculation 
of the MOKE can be downloaded from Ref. [H 

The small complex Kerr angle is generally defined as 
the off-diagonal divided by the diagonal Fresnel reflection 
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coefficients, for instance for incident p-polarized light: 
$p = r ps/ r pp — + ie. We would like to point out that 
this definition is insufficient as it does by no means define 
the sign of the Kerr angle, which depends on the choice 
of coordinate system, relative orientation of the p and 
s directions, the in-plane wave-vector of the light, and 
the sign of the exponent of the wave function. The sign 
convention in magneto-optics is indeed a long-standing 
problem, i.e. different authors report different signs for 
the complex Kerr and Faraday rotation angles^. Even 
worse hardly any article gives a clear definition of the em- 
ployed sign conventions. Therefore, we will here briefly 
derive an unambiguous definition of the sign, which must 
include the rotational sense of the Kerr angle, the geom- 
etry of the experimental setup, and the orientation of 
the longitudinal and transversal components mi and nit 
of the magnetization. We define the complex rotation 
to be positive when the rotational vector is pointing in 
the propagation direction of the reflected light, i.e. the 
polarization vector (p direction) is rotated in clockwise 
direction when looking in the direction of the reflected 
beam. The orientation of the coordinate system is de- 
fined as depicted in Fig. [1] with positive mi, which is also 
the direction of a positive external field pointing to the 
right, mt pointing up, and k\\, which is the direction of 
the in-plane light wavevector, pointing to the left when 
looking onto the the sample. With this convention for 
a 60 nm Fe film the Kerr rotation due to a positive m; 
results in a negative Kerr rotation 9 and a positive Kerr 
cllipticity e. The negative Kerr rotation is equivalent in 
sign to a counterclockwise turn (looking in direction of 
the incident beam) of the polarizer and the positive el- 
lipticity e to a clockwise turn of the polarizer out of the 
p direction, respectively. 

It is convenient to expand the Kerr angle <I> as a func- 
tion of the directional cosines of the magnetization vector 
M ^ e.g. for in-plane magnetization: 

$ = ^hmj+qi,i'm,i^rn t .i+q2^ml i +t i m ttl +0(m 3 ) , 

layers i 

(3) 

where U are the longitudinal (LMOKE), qij and q^^ the 
quadratic, and t, the usually much smaller transversal 
(TMOKE) response coefficients. Similar relations hold 
for the more general case including out-of-plane magne- 
tization, the Faraday effect, and even the calculation of 
Brillouin light scattering intensities^ which are closely 
related to the MOKE problem. 

Note that some authors use a different form for the 
second quadratic term, namely q2 (mf — m 2 ) , which is 
equivalent to our term, if the magnetization is in a single 
domain state except for a factor of —2 and a constant off- 
set of 1. In the case of a multi-domain state this kind of 
description does not generally hold as additional signifi- 
cant MO effects due to the magnetization gradient^ can 
be present. The longitudinal coefficients U stem from the 
linear MO coupling parameter K alone and are known to 
be isotropic, i.e. independent on the sample orientation, 



as long as the FM layers have cubic symmetry. By con- 
trast, the quadratic coefficients are due to a combined 
effect of the linear and quadratic MO couplings and are 
anisotropic, i.e. they depend on the relative orientation 
of the sample with respect to the plane of incidence. For 
cubic systems the resulting q coefficients have been found 
to have the form^ 

Qi = Qooi + (qon - 9ooi ) sin 2 (27) (4) 
32 = 2^011 - 3ooi ) sin(47), (5) 

where 7 is the angle between the in-plane component of 
the light wavevector and an in-plane Fe[001] direction, 
and gooi and go 11 are QMOKE constants for the plane 
of incidence parallel to the [001] and [Oil] directions, re- 
spectively. The anisotropy of the QMOKE stems from 
the symmetry of the effective SO coupling tensor Gijks, 
which is closely related to the symmetry of the crystal. 
For Fe(001) with in-plane magnetization (Gn — G12) and 
2G44 are the relevant SO coupling parameters for the 
[011] (7 = 45°) and [001] direction (7 = 0°), respec- 
tively. The parameter AG = (Gn — G\%) — 2G44 is a 
measure for the anisotropy strength^ 

The q\ coefficients have an isotropic and an anisotropic 
contribution. They depend on the sum of —K 2 /n 2 and 
Gn - G12 — AG cos 2 (27) 3 Therefore, errors in the de- 
termination of n and K affect the accuracy of the G val- 
ues. A wrong sign of the q\ value will not simply lead 
to wrong signs, but to wrong values of the second order 
MO coupling constants. On the other hand, the qi coef- 
ficients stem purely from the G parameters, namely from 
AG/2sin(47)ji£ and vanish if the plane of incidence is 
parallel to the symmetry directions [001] and [011]. 

The response coefficients can be determined experi- 
mentally - at least for single layer systems - with a suit- 
able setup, for instance, by rotating the fiekh 13 i 29 They 
can also be calculated numerically from the optical and 
MO material parameters n, K, and G. As our MOKE 
setup does not allow for a field rotation, we choose here 
the alternative approach of analyzing remagnetization 
loops recorded at different sample orientations. 

Typical MOKE loops for a 60 nm Fc film are shown in 
Fig. [2] The experimental Kerr rotation (top graphs) and 
cllipticity (bottom graphs) are plotted with connected 
open circles. We have determined the MOKE response 
coefficients I, q\, and q<i by fitting experimental remag- 
netization loops to a single domain model taking into 
account the sample orientation 7, the cubic anisotropy 
parameter K c /M s , and describing the Kerr angle via Eqs. 
([3H5]). While the left loops recorded at an angle 7 s» 45°, 
i.e. with field parallel to a hard [011] direction, depend 
on I and gon, the right loops are recorded at 7 rs 22.5° 
and, therefore, depend on I and both go 11 an d gooi- Thus, 
a simultaneous fitting of the loops for both orientations 
yields a full set of MOKE response coefficients I, gooi, 
and g ii. 
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hand, the imaginary part of the wavevector k± leads to 
a decreasing intensity of the electromagnetic wave with 
increasing depth inside the film, which determines the in- 
formation depth of about 40 nm, where the slope begins 
to asymptotically flatten. 



-0.1 0.1 -0.1 0.1 

Magnetic field (T) 

FIG. 2: Field dependence of measured (open circles) and cal- 
culated (red lines) Kerr ellipticity (bottom graphs) and Kerr 
rotation (top graphs) for a 60 nm Fe film. The left curves 
correspond to a sample orientation with hard axis parallel to 
external field (7 sa 45°), while for the right curves the field is 
parallel to an intermediate direction (7 ~ 22.5°). The field 
is swept form positive to negative values. The magnetization 
alignments are show with short arrows, the easy axis direc- 
tions of the magnetocrystalline anisotropy with crossed long 
arrows. 



IV. RESULTS AND DISCUSSION 
A. Linear MOKE 

The thickness dependence of the experimental linear 
MOKE extracted from the hard axis loops is marked in 
Fig. [3] by red circles and blue triangles for sample A and 
B, respectively. The data for thicknesses of 5nm and 
8nm measured for both samples are in excellent agree- 
ment indicating that the sample quality and thickness 
calibration of samples A and B are very similar. The up- 
per and lower curves corresponds to the imaginary and 
real part of the Kerr angle (ellipticity e and rotation 6) , 
respectively. In agreement with the results of previous 
publications on similar systems,— the Kerr ellipticity in- 
creases linearly with the film thickness below 5 nm, which 
is indicative for the law of additivity of the MOKE effect 
size to be valid in this regime. On the other hand, there 
is no such linear behavior for the Kerr rotation, which 
changes sign at about 4nm, meaning that the additiv- 
ity does not hold for the Kerr rotation. This breakdown 
of the additivity law is due to the dominant imaginary 
phase of the Kerr angle found for ultrathin Fe layers on 
noble metal substrates. In contrast, on semiconducting 
substrates, e.g. GaAs, the Kerr angle of Fe is mainly 
real and the additivity holds for the rotation, but not 
for the ellipticity^ For thicker layers, the phase of the 
electromagnetic wave inside the Fe layer changes due to 
the real part of the perpendicular wavevector component 
k± and therefore, gives rise to negative Kerr ellipticity 
contributions coming from Fe layers buried deeper than 
about 12 nm, where the slope changes sign. On the other 




20 30 40 50 

Fe thickness (nm) 

FIG. 3: Thickness dependence of the measured (red circles 
and blue triangles for sample A and B, respectively) and cal- 
culated (lines) linear ellipticity e (upper curves) and rotation 
6 (lower curves). Error bars are comparable to the symbol 
size. Black curves result from a calculation with material pa- 
rameters from literature. A fit to the data of sample A with 
thickness range - 8nm (red curves) yields a clearly differ- 
ent behavior than a fit to the data of sample B with a larger 
thickness range 5 - 60 nm (blue curves). Inset: Magnification 
of the low thickness region. 

The thickness dependence calculated using literature 
values of the indices of refraction from Refs. and 
El n Ag = 0.27 + 4.66i, n Au = 0.10 + 3.65i, n Fe = 
3.57 + 4.02i, and the linear MO coupling from Ref. HH 
Q = 0.0437 + 0.0040i is plotted as black lines. As found 
earlier by Qiu et al.^ the material constants from litera- 
ture approximately reproduce the Kerr ellipticity, which 
is insensitive to small phase changes of Q. However, the 
literature data fail to describe the Kerr rotation. 

We have fitted our experimental data employing the 
full 4x4 matrix formalis m 23 i 24 using fixed indices of re- 
fraction for Ag and Au from literature as specified above, 
and treating the index of refraction and the MO coupling 
of Fe as free parameters. The red and blue curves in Fig. 
[3]correspond to the data of samples A and B, respectively. 
The fit results are listed in Table |U It turns out that a 
satisfactory fit over the whole thickness range with only 
one thickness-independent set of material parameters is 
impossible. The fit to the data of sample A with smaller 
thicknesses results in a significantly about 10% larger in- 
dex of refraction with different phase as compared to the 
thicker sample B. On the other hand, the MO coupling 
parameter Q mainly differs in phase by about 10°, while 
the absolute values are in rather good agreement within 
less than 3% difference. Thus, it seems that the differ- 
ence between thin and thick Fe layers is mainly of optic 
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rather than of magneto-optic origin. 

Our value of the index of refraction of Fe determined 
from the data of the thicker sample B, n,p e = 3.53 + 3.72i, 
compares reasonably well with the value n = 3.57 + 4. 02i 
of Yolken and Kruger^ . While the real part is in excel- 
lent agreement, our imaginary part is about 7% smaller, 
which is probably within in the range of the systematic 
experimental errors. On the other hand, our value for 
the linear MO coupling Q for sample B is significantly by 
about 20% smaller and has a phase difference of about 
17° compared to the data of Krinchik and Artemevi^ 

The curve fitted to the data from the thinner sam- 
ple departs from the experimental data at about 15 nm, 
which corresponds to approximately half the penetration 
depth of the light. This circumstance might hint at an 
improper description of the optical properties of the Ag 
substrate as a reason for the apparent thickness depen- 
dence of the index of refraction of Fe. The substrate 
mainly influences the Kerr angle for Fe thicknesses be- 
low 15 nm as the light reflected from the substrate can 
reach the sample surface. In order to test this conjecture, 
we have additionally fitted the data with the indexes of 
refraction of the substrate and capping layers as free pa- 
rameters. However, we could not substantially improve 
the overall quantitative agreement of the fits. There- 
fore, improper optical parameters of the non-magnetic 
layers can be ruled out as a reason for the encountered 
thickness dependence of the optical Fe properties. Pos- 
sible explanations for the thickness dependence are: (i) 
The tensile strain of the Fe due to the small lattice mis- 
match of 0.7% between Fe and the Ag substrate leading 
to an anisotropic permittivity tensor, (ii) altered elec- 
tronic properties of the thin Fe layer due to the proxim- 
ity to the noble metal substrate and the capping layers, 
which can have a sizable influence^^ and (iii) interfa- 
cial MOKE contributions^ which have been neglected 
in the theoretical description. 



B. Quadratic MOKE 

The thickness dependence of the QMOKE coefficients 
qi for the hard axis (Fe[011] direction) and the easy axis 
(Fe[001] direction) configurations are plotted in Figs. 0] 
and 03 respectively. In contrast to LMOKE, in this case 
the real part (Kerr rotation 9) is the quantity, which 
depends linearly on thickness for ultrathin Fe layers, 
while the slope of the imaginary part (Kerr ellipticity e) 
changes sign at approximately 7 and 2nm for the [Oil] 
and [001] configuration, respectively. Thus, taking into 
account the significant QMOKE contribution the addi- 
tivity holds for neither phase of the Kerr angle. 

The extracted relevant second order MO coupling con- 
stants, Gn — G\2 and G44, for the hard and easy axis 
configurations are listed in Table |TJ They have been de- 
termined by fitting the data taking into account the n and 
Q parameters from the fits to the LMOKE data (Fig. [3]). 
The determination of the gooi parameter depends on the 
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FIG. 4: Thickness dependence of the measured (red circles 
and blue triangles for sample A and B, respectively) and cal- 
culated (lines) QMOKE for the Fe[011] direction (hard axis) 
parallel to the plane of incidence. 
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FIG. 5: Thickness dependence of the measured (red circles 
and blue triangles for sample A and B, respectively) and cal- 
culated (line) QMOKE for the Fe[001] direction (easy axis) 
parallel to the plane of incidence. 



both the values of gon and / [see Eq. ©]. Therefore, 
the propagation of errors results in a significantly poorer 
quality of the experimental data. As a consequence we 
could not reliably determine the G44 parameter for sam- 
ple A. 

We find maximum absolute values |gi| of the QMOKE 
at about 22 nm of 0.23 mrad and 0.37 mrad for the hard 
and easy axis configurations, respectively. These val- 
ues should be compared to the record QMOKE value 
of about 1.05 mrad^ recently found in Co2FeSi alloys 
in Ref. [H. This comparable order of magnitude of the 
QMOKE of Fe and Co2FeSi implies that the maximum 
QMOKE value of Fe in the visible wavelength region 
might be even larger than that of Co2FeSi as both materi- 
als are expected to have a distinct frequency dependence 
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TABLE I: Index of refraction (n), linear (Q) and quadratic (Gn — G12, G44) MO coupling constants derived from data of 
samples A and B. 



Parameter 


Sample A (0 - 8 nm) 


Sample B (5 - 60 nm) 


Literature 2 ^ 


n 


4.06 ±0.03 + (3.85 ± 0.03)i 


3.53 ±0.03 + (3.72 ±0.03)i 


3.57±4.02i 


Q 


0.0331 ± 0.0002 - (0.0127 ± 0.0002)i 


0.0356 ± 0.0004 - (0.0074 ± 0.0003)i 


0.0437 + 0.004i 


Gn — G12 


-0.0544 ± 0.0005 - (0.0287 ± 0.0005)i 


-0.0358 ± 0.0002 - (0.0382 ± 0.0002)i 




G44 




-0.0117 ± 0.0003 - (0.0349 ± 0.0003)i 





resulting from their electronic band structures. 

In Fig. 2 of Ref.|l3| Postava et al. give the dependence 
of LMOKE and QMOKE contributions on sample orien- 
tation for a 50 nm thick bcc-Fe(OOl) sample capped with 
1.5 nm Pd and measured at an incident angle of 3.25°. 
Based on our fitted optical and MO material constants 
we have calculated the Re(l), i?e(qooi)i & n d i?e(qon) con- 
stants for the sample structure and experimental config- 
uration of Ref.[l3| and assuming npd = 1.87 + 4.44i. We 
find a value of Re(l) = — 0.066 mrad which is compara- 
ble, but about 20% smaller than the value of Re (I) ~ 
—0.09 mrad— determined by Postava et al. On the other 
hand, our QMOKE data differ more distinctly, although 
they are of a similar order of magnitude. While we find 
Re{qon) = 0.079, which is about 30% larger than the 
value Re(qon) ~ 0.06 determined by by Postava et al, 
our value for -Re(gooi) = —0.073 is about a factor of 2 
smaller than J?e(gooi) ~ —0.16 determined by Postava 
et al. This difference of the QMOKE is clearly larger 
than the experimental errors caused by uncertainties in 
the layers thicknesses, the optical properties of the cap- 
ping layers, and the calibration of the MOKE setup. Al- 
though both the sample of Postava et al. and our samples 
are epitaxial bcc-Fe(OOl), a possible explanation could 
be a strong structural dependence of the QMOKE as de- 
scribed in Ref. for Co 2 FeSi. 



V. CONCLUSIONS 

The thickness dependence of the linear and quadratic 
MOKE of wedge-type thin Fe(001) films magnetized in 



the sample plane has been measured. Good quantitative 
agreement of the experimental data with calculations as- 
suming bulk-type optic and MO material constants in- 
dicates that the thickness dependence of the MOKE is 
mainly due to bulk-type magnetooptical coupling. On 
the other hand, we found a sizable departure from theory 
for thicknesses below about 10 nm. This can be explained 
by thickness dependent optic and MO material param- 
eters, which are possibly due to MO surface effects or 
thickness dependent features of the bandstructure, e.g. 
quantum-well states. By fitting the experimental data 
to results of a numerical model, we extracted a complete 
set of material constants n, K, (Gn — G12), and G44 for 
the quantitative description of the MOKE of bcc-Fe(OOl) 
at the laser frequency employed. To our knowledge this 
is the first report of the second order MO-coupling con- 
stants of Fe. They are comparable to the first-order con- 
stants and thus, of general significance for the theoretical 
description of the MOKE. The index of refraction n is 
in excellent agreement and the linear MO-coupling con- 
stants K agree reasonably with previous works. In con- 
trast, a comparison of the second-order constants with 
earlier QMOKE data from Postava et al. demonstrates a 
remarkable difference. This suggests that the anisotropic 
second-order MO-coupling might strongly depend on the 
sample properties. 
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